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Jo ■ Abstract 



We investigate the FIR-absorption of short period parallel quantum wires 
in a perpendicular quantizing magnetic field. The external time-dependent 
electric field is linearly polarized along the wire modulation. The mutual 
Coulomb interaction of the electrons is treated self-consistently in the ground 
q ■ state and in the absorption calculation within the Hartree approximation. We 

consider the effects of a metal gate grating coupler, with the same or with 
a different period as the wire modulation, on the absorption. The evolution 
of the magnetoplasmon in the nonlocal region where it is split into several 
Bernstein modes is discussed in the transition from: narrow to broad wires, 
^) ■ and isolated to overlapping wires. We show that in the case of narrow and 

not strongly modulated wires the absorption can be directly correlated with 
the underlying electronic bandstructure. 
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I. INTRODUCTION 



Parallel quantum wires can be produced in ALGai-aAs-GaAs heterostructures by a 
periodically modulated metal gate on top of the Al^Gai^As insulator. By applying a 
voltage bias between the gate and the two-dimensional electron gas (2DEG) at the insulator- 
semiconductor interface the properties of a 2DEG can be studied in the entire range from a 
homogeneous system to isolated wires. 

The theory of far-infrared (FIR) absorption of the homogeneous 2DEG in the presence 
of the metal gate acting as a grating coupler has been developed by Zheng, Schaich, and 
MacDonald! and by Liu and Das Sarmap. In those papers the FIR response has been 
calculated semiclassically, allowing the presence of a weak magnetic field. The FIR absorp- 
tion of a modulated system, for instance due to the presence of a gate potential, has been 
extensively studied experimentally for a simple unidirectional modulation! and calculated 
self-consistently for a sinusoidal potential and strong magnetic fields.!'! Modes specific to 
more structured modulation like parallel pairs of wires have been identified as well.! 

The well known magnetoplasmon dispersion law for a homogeneous 2DEG in a strong 
magnetic field is 

to = lu c + 2nn s e q/ (nm) , (1.1) 

which is valid in the lowest order in the wave vector q, i. e. as long as the second term 
is much smaller than the first one. For shorter wavelengths or for lower magnetic fields 
the magnetoplasma oscillations have higher branches, around harmonics of the cyclotron 
frequency, noo C} n = 2,3, known as Bernstein modes.! The Bernstein modes are likewise 
known in single wires and dots.!'! 

In this publication we want to focus our attention on the FIR absorption of a periodically 
modulated 2DEG in a perpendicular quantizing magnetic field, in the regime where the 
Bernstein modes are a prominent feature of the frequency dispersion. Furthermore, we take 
into account the grating coupler effects of the metal gate on the incident FIR radiation. 
The external time-dependent electric field is hence modulated with a large wavevector, the 
lowest inverse lattice vector of the gate structure. This is an approximation to the leading 
order compared to the studies of the grating coupler effects mentioned above.!'! In principle, 
the 2DEG can be modulated by an independent method, i. e. not by the gate structure. 

We calculate the FIR absorption fully quantum mechanically with the help of the inverse 
dielectric function. The mutual Coulomb interaction of the electrons is treated within the 
Hartree approximation both in the ground state and in the absorption calculation, i. e. the 
absorption is obtained in the framework of the random phase approximation (RPA). The 
intra- or inter-wire Coulomb interactions thus enter the model on equal footing. 

We show that as the width of the wires is reduced the resulting Landau bandstructure can 
be identified through the effects of its van Hove singularities (vHS) on the FIR absorption. 
The evolution of the FIR absorption is traced from the case of weakly coupled wires to 
the case of a weakly modulated 2DEG by increasing the electron density but keeping the 
modulation strength constant. We systematically use the induced density to sort out the 
complex hierarchy of absorption modes caused by the modulation of the electron gas or by 
the electric coupling between wires. 



2 



II. MODEL 



We consider a modulated 2DEG located in the plane r = (x, y), and describe it with the 
Hamiltonian 

H = H + V mod + V H . (2.1) 

H Q = (P+eA) 2 /2m is the Hamiltonian of the noninteracting system. We employ the Landau 
gauge for the vector potential, A = (0,Bx). K tlo d(x) = V cos(Kx) is the simplest model 
of an electrostatic potential modulation varying in only one direction. Vh is the Hartree 
potential felt by each electron, self-consistently with the particle density, 

P( x ) =J2Pp cos(pKx) , (2.2) 



p>0 



e f j_/ P( x> ) 



V H (x) = —[ dr 

K J 

e 2 2tt ^ p p 



r — r 



V^oosCpifx), (2.3) 
p>i P 

where k is the dielectric constant of the semiconductor host of the 2DEG. We assume 
the electrical neutrality to be ensured by a uniform background of positive charges, of 



density — po- The eigenstates of the Hamiltonian (|2.1|) are calculated numerically, within 
an iterative scheme, by successive diagonalization in the basis of the wave functions cor- 
responding to Hq, known as Landau wave functions. The resulting single particle ener- 
gies, E nt x = E nt x +a, a, = 2ir / K , have a one- dimensional periodic band structure, where 
n = 0,1, 2, ... is the Landau quantum number and Xq is the center coordinate. We assume 
spin degeneracy. The eigenfunctions of H have the form L~ 1//2 exp(— iX y / l 2 )ip n:Xo (x) , where 
L is a normalization length and I = (h/eB) 1 / 2 is the magnetic length. 

We consider a time- dependent electric field incident to the 2DEG, linearly polarized 
in the direction of the modulation, x, with only one Fourier component, of wave vector 
q = (q, 0) and frequency u. This field is supposed to be sufficiently small, within the linear 
response regime. We denote by (f> eX t(q,uj) the associated electric potential. 

In order to evaluate the absorption due to plasma oscillations we need to calculate first 
the dielectric matrix £gg'(q, u), where G and G' are vectors in the reciprocal space, of the 
form mK, m = 0,±1,±2, .... We use the random-phase approximation, as described e.g. 
in the paper by Wulf et. al.,@ in which e GG "(g, cu) = 5gg' ~ -^^\ Xgg'{q, uj )- The dielectric 
susceptibility Xgc is given by the Lindhard formula, which in our case reads 

1 ^ r a 



nn' 

) - T{E n , 

,X ) 



E n ,Xo - Eri,x -Tko -ir) 

X Jnn'XoiQ + G)JL';xM + , (2.4) 

where T is the Fermi function, 77 — > + , and Jnn'Xo^) = (Vv Xol^^l^PnXo) ■ 

The absorbtion power can be calculated from the Joule law of heating, which may be 
written as@ 
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xg|0 cxt (g,^)| 2 , (2.5) 

where q = qi + G. Due to the periodicity of the system SQQ'(q, uj) = £g-k,g'-k(<1 + K,uj), 
and we can take < q± < K. We will be mostly interested in the case when the wavelength 
of the incident field is identical to that of the modulation, as in the absorption experiments 
with grating-coupler devices, that is q = q\ = K and G = 0. For simplicity, we will normalize 
the external potential in Eq. Q2.5| ) such that ^g|0 ext (g, cu) | 2 = 1. 

In order to avoid calculation of the plasma poles of the dielectric matrix, we assume 
a finite i] = hu c /50, in Eq. ( |2.4j ). In other words we assume an unspecified dissipation 
mechanism, which leads us directly to the oscillator strength of the collective modes, and 
thus to measurable results. 



III. RESULTS 

The parameters of our model are those for GaAs, the electron mass m = 0.067m and 
k = 12.4. The temperature will be fixed to 1 K. We truncate the Fourier series such that 
| mK \< MK, and we achieve numerical convergence for 5 < M < 15. 

As has been mentioned in the introduction the magnetoplasmon dispersion at low mag- 
netic fields and short wavelengths splits into Bernstein modes due to the interaction with 
harmonics of the cyclotron resonance 

In Fig. la we show the absorption power in the regime of the Bernstein modes, for 
the homogeneous system, for wave vectors varying from O.Oli^ to K, with K = 27c/a, and 
a = 50 nm. The electron density is fixed to po = 1.8 x 10 11 cm" 2 and the magnetic field is 
B = 3.74 T, which corresponds to the filling factor v = 2. Eq. (|1 . 1| ) holds only for wave 
vectors below q = 0.2K for the present choice of parameters. In Fig. lb we display the 
results for a modulated system, with a modulation period a = 50 nm, and a modulation 
amplitude V = 5 meV. As we observe, each peak corresponding to a Bernstein mode splits 
into several subpeaks. The same modulation effect is evidenced in Fig. lc, where we have 
fixed the wave vector of the incident field, q = K, but we sweep the magnetic field such that 
the filling factor varies between 1 and 3. The traces for u = 1.2 and v = 2 are magnified in 
Fig. 2. Each Bernstein mode has up to four subpeaks. 

We start to analyze this complicated internal structure, by showing in Fig. 3 the energy 
spectra corresponding to Fig. 2. The periodic Landau bands are displayed in half of a 
Brillouin zone, i.e. for < KXq < it. Due to the short period, in the absence of the 
Coulomb interaction the energy bands are not parallel. This is seen here only for v = 2, 
where the Fermi level is in an energy gap, and thus the screening effects are small and the 
energy dispersion large, Fig. 3a. For v = 1.2, Fig. 3b, the Hartree screening reduces the 
energy dispersion to very narrow bands. 

The absorption peaks of Fig. 2 reflect the nonparallel Landau bands, or, equivalently, 
a certain dispersion, in the center-coordinate space, of the energy interval between pairs 
of Landau bands. This is best seen for v = 2, where for each Bernstein mode we can 
energetically relate the lowest and the highest subpeaks to the inter-Landau-level transitions 
around the center and around the edges of the Brillouin zone, i. e. X = and X = n/K 
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respectively. The energy dispersion is locally flat in those regions, or in other words the 
density of states has van Hove singularities, such that the single-particle transitions may 
turn into collective excitations with a slightly different energy. 

For instance, for 2.%uj c < uj < 3.6u c , the peak at uj = 2.95u c corresponds to the energy 
interval £3,0 — -Eb,o — 2.77uj c , plus a blue shift determined in part by the increased energy 
distance between the adjacent levels, E 3 ^x — E 0) x , with say < X < 0.5/ K, and in part 
by the external electric field, similarly to the second term of Eq. ( |1 . 1| ) . The highest peaks, 
at uj = 3.32co> c and uj = 3.37u; c , apart of this small splitting, can be related to the energy 
interval between the other pair of vHS, E 3 ^/k — Eq^ik = 3.27uj c . The small splitting of 
the highest mode cannot be explained only in terms of the energy spectrum. The middle 
peak, at uj = 3.18u; c , may be put in correspondence with the average energy interval, and 
thus considered as a combination of the two types of vHS modes. Qualitatively, the same 
structure is obtained for the Bernstein modes around the other multiples of uj c . 

For v = 1.2, the dashed line of Fig. 2, the states around the central vHS of the Brillouin 
zone become empty, and the corresponding subpeaks vanish. The bands are nearly flat, due 
to the screening effect, but the splitting of the higher subpeaks is still present. Similar results 
have been very recently obtained by Brataas et. alH for the singularities of the imaginary 
part of the dielectric function of modulated systems, by analytical calculations, without the 
inclusion of the screening in the groundstate, but also related to the nonparallel Landau 
bands. 

In order to support the above interpretation of the Bernstein subpeaks, in Fig. 4 we 
show the particle density at equilibrium, p(x), the induced density, pi n d(x), and the extreme 
density profiles during the oscillations, p(x) ± Pi n d(x), inside a unit cell, < x < 2tt, for 
v = 2 and q = K. Obviously, the minima (maxima) of the equilibrium density correspond 
to the maxima (minima) of the modulation-potential energy. For graphical reasons we 
have amplified the induced density. The frequency of Fig. 4a, uj = 2.70uj c , is chosen 
between two collective modes, where the absorption power is very small. The characteristic 
feature of such absorptionless oscillations is that the evolution of the density preserves the 
reflection symmetry of the unit cell. A kind of a breathing mode, attributed to single-particle 
excitations in the presence of a dissipation mechanism. 

For the frequencies corresponding to collective peaks, Fig. 4b and c, we observe a global 
shift of the particle density inside the unit cell, i. e. we observe dipolar modes. In Fig. 4b both 
the density minima, around Kx = and Kx = 2ir, and the maxima, around Kx = n, have 
horizontal oscillations along the x direction. The amplitude for the minima is larger than 
for the maxima, which is in agreement with the domination of this mode by the transitions 
around the center of the Brillouin zone (associated to the maximum modulation-potential 
energy). We also observe a strong breathing component, due to the continuous spectrum of 
single-particle transitions of higher energies, which can be associated with the blue shift. In 
Fig. 4c we have dipolar oscillations of the density maxima, i. e. of frequency corresponding 
to the lateral vHS, with only a small breathing contribution, due to the excitation gap above. 
The neighboring mode, with uj = 3.32u; c , Fig. 2, is very similar to that of Fig. 4c, but having 
an additional node in the induced density. Also, the soft mode for uj = 3.18uj c is similar to 
that of Fig. 4b, with one node less in the induced density, but with a stronger breathing 
component. 

We will discuss now the situation when the modulation period is much larger than the 
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magnetic length. In this regime the Landau bands are parallel, and the screening effects are 
stronger than before. In Fig. 5 we display the absorption for frequencies within the interval 
0.5u; c < uj < 3u) c , again with a complicated structure of the Bernstein modes. As in the 
case with a = 50 nm. this structure also reflects the many branches of the magnetoplasma 
oscillation spectrump but now the internal peaks cannot be related to the energy spectrum 
as before, because the excitation energies between the central and the lateral vHS coincide, 
e. g. like in Fig. 6. Here the energy dispersion in the vicinity of the Fermi level is flat due 
to the electrostatic screening. Ill Since the states around the center of the Brillouin zone are 
less populated than those at the edges, the active frequencies in the absorption spectrum 
correspond to dipolar oscillations of the density maxima, as shown in Fig. 7. The absorbant 
modes differ by the number of nodes in the induced density. In particular, for v = 1.6, the 
height of the absorption peaks is lower for the modes with fewer nodes, which is somewhat 
contrary to the expectation. 

A similar situation is shown in Fig. 8. Here we have reduced the mean particle density to 
Po = 0.5 x 10 11 cm -2 , in order to explore the limit of isolated, parallel quantum wires. The 
magnetic field is fixed to B = 4.67 T. Even if the strips inbetween the wires are depleted to 
nearly zero particle density in the ground state, the wires still interact in the presence of the 
incident electric field. Consequently, the oscillations with a high number of nodes, capable of 
penetrating into the depleted strips, like in Fig. 8b, are energetically favorized with respect 
to those with less nodes, but which keep a rigid wire separation, Fig. 8c. Here we depict the 
density configurations for the middle and the right peaks in the group with three maxima 
of Fig. 8a. For strictly isolated wires the induced density would be either symmetric or 
antisymmetric with respect to the wire center. This condition is the best obeyed by the 
mode at uo = 1.26u c , which is thus specific to the wire response, the other modes being still 
under the influence of the electric coupling between the wires. 

In Fig. 9a we display the evolution of the absorption spectra with increasing density, 
from the limit of isolated wires, po = 0.1 x 10 11 cm -2 , up to well coupled wires, po = 0.9 x 10 11 
cm -2 . The wave vector of the time dependent electric field is now small, q = 0.1K, in order 
to concentrate the attention only on the effects of increasing the electron density on the 
absorption in a periodic system. The effect of the gate, that we neglect here, is mainly to 
repeat the same for each Bernstein mode. The Fermi level is always in the Landau band 
with n = 0. For the lowest density the energy dispersion is large, but only few states are 
occupied at the bottom of the energy spectrum, while for the highest density the Landau 
band becomes very narrow, due to screening, such that the depleted strips vanish. In the 
limit of isolated wires the induced density at the single absorption peak has only two nodes 
inside each unit cell, as can be seen in Fig. 9b. This kind of motion has thus a strong 
breathing component, i. e. vertical oscillations of the total charge of the type shown in Fig. 
4, and only small dipolar (horizontal) components which vanish with changing frequency. 
Hence the internal motion of the electrons inside the wires has a strong contribution to 
the resonant modes, or in other words Kohn's theorem is not obeyed here, due to the 
nonparabolic, soft lateral confinement. 

Increasing the density, the wire coupling allows additional charge oscillations inside the 
unit cell, Fig. 9c, which may cumulate a stronger total dipolar component, yielding thus 
higher absorption peaks. The lower the frequency, the higher the number of nodes within 
each wire, as also noticed for q = K, in Fig. 8. The mode at u = 1.26o; c , with few nodes per 



6 



cell, corresponds to the only mode of the isolated wire system. Such a mode is present in 
each trace of Fig. 9a, being the highest in frequency but the lowest in amplitude. Instead, 
the new modes, with additional internal oscillations, grow rapidly at lower frequencies, when 
the density, and thus the wire coupling, increases. 

For a given mode, the characteristic wavelength of the internal motions inside the unit cell 
is determined by the lateral confinement and by the wire coupling. In our periodic system 
both of them result self-consistently from the external modulation potential, the magnetic 
field, and the electron density. The evolution from a single-peak absorption around u c , spe- 
cific to isolated wires, to the multi-peak structure specific to the modulated system, occurs 
thus by the progressive appearance of new modes, with richer internal motion, and accu- 
mulating larger oscillator strength. For each density, the number of nodes in the absorption 
maxima of Fig. 9a gradually increases from right to left. When the density is increased 
further, beyond what is shown in Fig. 9, such that more than one Landau band is occupied 
a stable single peak emerges with energy and induced density specific to the homogeneous 
system. 



IV. SUMMARY 

In the present publication we have explored the properties of the FIR-absorption of a 
2DEG modulated in one direction, varying the modulation and the electron density in order 
to cover the whole range from isolated but electrically coupled wires to a weakly modulated 
2DEG. The electron-electron interaction has been treated self-consistently within the time- 
dependent Hartree approximation. Experimentally the transition from wires to an almost 
homogeneous 2DEG can be achieved in heterostructures with a modulated metal gate, which 
in turn also acts like a grating coupler for the incident FIR radiation. The modulation of 
the incident field should be determined self-consistently including effects from the gate and 
the 3D structure of the sample.^ To avoid this problem that is further complicated by 
the presence of the magnetic field we considered primarily the effects of the two lowest 
Fourier components in the incident field, (q « 0, q = K). The presence of the short 
wavelength component and the magnetic field of intermediate strength place the dispersion 
of the magnetoplasmon in the regime of Bernstein modes. 

Our calculations show that information about the nonparallel Landau bands in short 
period wires can be extracted from the peak structure of the FIR absorption. 

The effects of the vHS of the underlying energy band structure on the dielectric function 
in weakly magnetically and electrically modulated system have been studied in absence of 
the Coulomb interaction between electrons.0'0 
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FIGURES 




FIG. 1. (a) Absorption spectra for a homogeneous system. K = 2ir/a, a = 50 nm. The filling 
factor is 2 and the density po = 1.8 x 10 11 cm~ 2 (b) The same as (a), with a modulation of period 
a = 50 nm and amplitude V = 5 meV. (c) The same modulation, with a variable magnetic field 
such that the filling factor evolves between 1 and 3. q = K. 
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FIG. 2. The curves of Fig. lc corresponding to v = 2, full line, and z/ = 1.2, dashed line, 
magnified. 
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FIG. 3. Energy spectra for (a) v = 2 and (b) z/ = 1.2 in half of the Brillouin zone. The dashed 
horizontal lines show the Fermi level. 
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FIG. 4. Density profiles for various indicated frequencies taken from the absorption spectrum 

of Fig. 2, v = 2. The full lines show the groundstate density, the dashed lines the maximal induced 
density, and the marked curves show the extreme density configurations. 




FIG. 5. Absorption spectra for a modulated system with a = 200 nm and V = 20 meV. The 
average density is po = 1-8 x 10 11 cm~ 2 , q = K. 
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FIG. 6. Landau bands for the modulation of Fig. 5, with v = 1.6. 
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FIG. 8. (a) Absorption peaks for the same modulation as for Fig. 5, but with the average 
density po = 0.5 x 10 11 cm -2 , (b) and (c) two density diagrams. 
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FIG. 9. (a) Absorption spectra for various electron densities, from 0.1 to 0.9 units of 10 11 cm~ 2 , 
covering the transition from isolated wires to a modulated system. The magnetic field is B = 4.67 
T and the wave vector of the incident electric field is q = 0.1K. (b) and (c) The induced electron 
density, in a full period of the perturbed system, for po = 0.1 x 10 11 cm~ 2 and po = 0.5 x 10 11 
cm~ 2 , for the frequencies corresponding to the absorption maxima. 
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